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ABSTRACT 

The possibility that the masses M, of supermassive black holes (SBHs) correlate with the total gravitational 
mass of their host galaxy, or the mass Mqm of the dark matter halo in which they presumably formed, is inves- 
tigated using a sample of 16 spiral and 20 elliptical galaxies. The bulge velocity dispersion cjc, typically defined 
within an aperture of size R <^ 0.5 kpc, is found to correlate tightly with the galaxy's circular velocity Vc, the latter 
measured at distances from the galactic center at which the rotation curve is flat, /? ^ 20-80 kpc. By using the 
well known M, - CTc relation for SBHs, and a prescription to relate Vc to the mass of the dark matter halo Mdm 
in a standard ACDM cosmology, the correlation between (Jc and v^. is equivalent to one between M, and Mom- 
Such a correlation is found to be nonlinear, with the ratio M./Mdm decreasing from 2 x lO"'* for Mom ^ 10'"^ M©, 
to 10"^ for MjyM ^ 10'^ Mq. Preliminary evidence suggests that halos of mass smaller than ~ 5 x 10" Mq are 
increasingly less efficient - perhaps unable - at forming SBHs. 



1. INTRODUCTION 

Never has the study of supermassive black holes (SBHs) been 
greeted with such attentiveness as in the past few years. Al- 
though by the mid 1960s it was generally recognized that the 
energy source of the newly discovered quasars must be gravi- 
tational (e.g. Robinson 1965), for the following three decades 
the existence of SBHs was destined to be surrounded by skep- 
ticism. By the mid 1990s, a few years after the launch of the 
Hubble Space Telescope, it was widely accepted (Kormendy & 
Richstone 1995). Today, SBHs feature prominently in models 
for the formation and evolution of their host galaxies (e.g. Silk 
& Rees 1998). Freed from the burden of having to demonstrate 
the very existence of SBHs, we can now begin asking cardi- 
nal questions: how did SBHs form, how do they accrete, how 
do they evolve, and what role do they play in the formation of 
cosmic structure? 

Tackling such questions has been made possible by the dis- 
covery that SBHs are tightly related to the large-scale properties 
of their host galaxies. Black holes masses, M,, correlate with 
the blue luminosity Lb of the host bulge (Kormendy & Rich- 
stone 1995). Even more compelling is the correlation between 
M, and the bulge velocity dispersion, M, oc cr" (Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). Using the most up-to- 
date set of SBH mass measurements such a correlation can be 
written as: 



M. =(1.66 ±0.32) X 10Vc/200kms"') 



K4.58±0.52 



Mr. 



(1) 



(Ferrarese 2002), where ctc is normalized to an aperture of size 
one eighth the bulge effective radius. The M, - CTc relation has 
supplanted the M,-Lb relation as the tool of choice for the 
study of SBH demographics (Ferrarese 2002; Yu & Tremaine 
2002), for the following reasons. The M,-Lb relation suffers 
from a large scatter (a factor of several in M,, Ferrarese & Mer- 
ritt 2000), most (but not all) of which is ascribable to lenticular 
and spiral galaxies (McLure & Dunlop 2002). The M, - Cc rela- 
tion, on the other hand, appears to be intrinsically perfect (Fer- 
rarese & Merritt 2000; Merritt & Ferrarese 2001a; Ferrarese 
2002). It holds true for galaxies of disparate Hubble types (from 



SBs to compact elUpticals to cDs) belonging to wildly different 
environments (from rich clusters to the field), showing perfectly 
smooth (e.g. NGC 6251) or highly disturbed (e.g. NGC 5128) 
morphologies. Furthermore, cr is a more faithful gauge of mass 
than Lb'- the M, -Cc relation thus betrays the existence of a 
tight connection between the masses of SBHs and those of the 
hot stellar component surrounding them (Merritt & Ferrarese 
2001b). 

Have we indeed found, in the M,-ac relation, the fossil 
record that links the formation and evolution of SBHs to those 
of their host galaxies? The M, -CTc relation, and the equally 
tight relation recently discovered by Graham et al. (2001) be- 
tween M, and the central concentration of bulge light, probe 
the connection between SBHs and the baryonic component of 
galactic bulges. Measurements of a typically do not extend 
far enough to penetrate the region dominated by the dark mat- 
ter (DM) component (Faber & Gallagher 1979). And yet, it is 
not the mass of the bulge but rather, the total mass M,ot of the 
galaxy or the mass Mom of the DM halo, which is the key ingre- 
dient in most theoretical models proposed for the formation of 
SBHs (Adams, Graff & Richstone 2000; Monaco et al. 2000; 
Haehnelt, Natarajan & Rees 1998; Silk & Rees 1998; Haehnelt 
& Kauffmann 2000; Cattaneo, Haehnelt & Rees 1999; Loeb & 
Rasio 1994). Once a correlation with M,,,, {or: Mqm) is recov- 
ered by the models, the correlation with bulge mass is implicit 
because, in standard CDM scenarios, the bulge mass is loosely 
determined by the halo properties (e.g. van den Bosch 2000; 
Haehnelt, Natai-ajan & Rees 1998; Zhang & Wyse 2000). 

This paper pursues, from an observational standpoint, the 
possibility of a link between M,, M,o, and Mom- A tight cor- 
relation is indeed found in the form of a relation between M, 
and the circular velocity of the SBH host galaxy, measured well 
beyond the optical radius. 

2. DATABASE 

Measuring the total gravitational mass of a galaxy is a 
formidable challenge. In bright ellipticals, the temperature and 
density distribution of the thermal X-ray gas can lead to a mea- 
sure of M,oi under the minimal condition of hydrostatic equilib- 
rium (e.g. Fabricant, Rybicki & Gorenstein 1984). However, 
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while much progress is expected as Chandra data are collected, NGC 891, Rogstad & Shostak 1972; NGC 1808, Saika et al. 
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Fig. 1 . — Correlation between bulge velocity dispersion ac and disk circular 
velocity Vc for a sample of 37 galaxies with either optical (open circles) or 
HI (solid circles) rotation curves. Points marked by errorbars only correspond 
to galaxies for which the rotation curves does not extend beyond R25. The 
galaxy to the far left, with the smallest value of Uc, is NGC 598. The solid Une 
corresponds to a fit to all galaxies with ac > 70 km s"' and R(vc)/R25 > 10. 
The dotted line corresponds to the fit to all galaxies, with the only exclusion of 
NGC 598. 



previous generations of X-ray satelUtes lacked the necessary 
combination of spectral response and sensitivity to allow an ac- 
curate measurement of all required observables (e.g. Biermann 
et al. 1989; Awaki et al. 1994; Matsushita et al. 1998). Stars, 
globular cluster and planetary nebulae (PNe) are good tracers 
of the gravitational potential (e.g. SagUa et al. 2000; Cote et al. 
2001 ; Amaboldi et al. 1998). However, the derived masses suf- 
fer large uncertainties due to the unknown velocity anisotropy 
of the system and, in particular for the case of PNe, the small 
sample size of the tracers (Napolitano et al. 2001; Merritt & 
Tremblay 1993). 

In spiral galaxies, however, the rotation curves of the cold 
disk component bear the direct imprint of the mass distribu- 
tion. Indeed, the lack of a decline in the circular velocity Vc 
at radii of order and beyond the optical radius of spirals was 
one of the first indications of the existence of a significant dark 
matter component (Rubin, Ford & Thonnard 1978, 1980; Ru- 
bin, Thonnard & Ford 1977; Krumm & Salpeter 1979). To first 
order, therefore, a relation between M, and M,o, would be re- 
flected in a relation between M, and Vc- 

Of the four spiral galaxies with a secure measurement of M, 
(see the compilation by Merritt & Ferrarese 2002) only two, 
the Milky Way and NGC 4258, have published rotation curves, 
both derived from 21-cm studies (OlUng & Merrifield 1998; 
van Albada 1980). However, the sample can be augmented if 
one is willing to exchange M, for the bulge velocity dispersion 
CTc, following equation (1). Of all spirals for which a has been 
measured, 46 have rotation curves, either in HI (21 objects) or 
optical (25 objects). For six of the galaxies with an HI rotation 
curve, warnings are found in the original papers as to the use of 
the published circular velocities for mass measurements. Rea- 
sons include warps in the HI disk preventing one from depro- 
jecting the observed rotation curve (NGC 4565, Rupen 1991; 



Fig. 2. — Residuals from the best fit to all galaxies with the exclusion of 
NGC 598 (dotted hne in Fig. 1), plotted as a function of R(vc)/R25. 



1980, Koribalski et al. 1993) and/or the presence of signifi- 
cant non-circular motion in the HI gas (NGC 3079, Irwin & 
Seaquist 1991; NGC 3031, Rots & Shane 1975; NGC 4736, 
van der Kruit 1974). For two of the galaxies with an optical 
rotation curve (NGC 3200 and NGC 7171), two values of a 
discrepant at the > 99.9% level exist in the literature (McEl- 
roy 1995; Whitmore et al. 1984; Dressier 1984); these galaxies 
have also been excluded from the sample. 

Data and references for the remaining 38 galaxies are given 
in Table 1. To allow for the application of equation (1), the 
bulge velocity dispersion a is corrected to the equivalent of an 
aperture of radius r,, / 8 following the prescriptions of Jorgensen, 
Franx & Kjaergaard (1995). The bulge effective radius is 
taken from Baggett et al. (1998, but see notes to Table 1). 

3. THE Vc-<Tc RELATION 

Fig. 1 shows the circular velocity versus the bulge velocity 
dispersion within rg/S for all galaxies in Table 1. The corre- 
lation is immediately apparent, but a few observations are in 
order 

First, a glance at Table 1 will reveal that the most noticeable 
difference between the two samples of optical and HI rotation 
curves is, not surprisingly, that while only two of the optical ro- 
tation curves extend beyond R25, the galaxy radius at the B = 25 
mag arcsec"^ isophote, all but one of the HI rotation curves 
do. Inspection of the original papers shows that the HI rota- 
tion curves all exhibit a common pattern. It is well known (e.g. 
Sofue & Rubin 2001 and references therein) that as one moves 
away from the center, rises rapidly, peaks somewhere be- 
tween 3 and 15 kpc, and then remains flat (sometimes after a 
modest initial decline) until the last measured point, typically 
Rivc) = 30-40 kpc. Optical rotation curves rarely extend be- 
yond 10 kpc; at this distance has not yet settled into a flat 
rotation curve. In Fig. 1, galaxies with rotation curves not ex- 
tending beyond R25 (shown as errorbars only) appear to dis- 
play larger scatter, relative to the mean correlation defined by 
aU galaxies, than the galaxies with a more extended rotation 
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Fig. 3. — Same as Fig. 1, except solid dots represent spiral galaxies with 
i?(Vc)//?25 > 1-0, open dots represent the sample of elliptical galaxies from 
Kronawitter et al. (2000). The dotted line is a fit to all spiral galaxies with 
the exception of NGC 598, while the solid line is the fit to all spirals with 
(Tc > 70 km s~' . The latter fit is consistent with the one obtained for the ellip- 
tical galaxies alone (see fits 3 and 4 in Table 2). The dashed line is a fit to the 
entire sample of ellipticals plus spirals with dc > 70 km s~' . 

Second, a simple visual inspection of Fig. 1 shows that 
among the galaxies with R(vc)/R25 > 1, NGC 598 (= M33, the 
galaxy with the lowest value of CTc in the plot) is a clear out- 
lier. Furthermore, the two galaxies with <Tc < 70 km s"' seem 
to have a significantly larger circular velocity that would be in- 
ferred if a Unear fit (in logCTe - log Vc) were to be performed to 
the rest of the sample. 

To quantify the above statements, we performed a regression 
analysis (Akritas & Bershady 1996) to all galaxies except NGC 
598 (Table 2 and dotted Une in Fig. 1). Fig. 2 shows the resid- 
uals from the best fit as a function of /?(Ve)//?25, confirming the 
larger scatter for galaxies with rotation curves not extending 
beyond .^25- For the scope of this paper, the usefulness of Vc is 
directly related to how confidently Vc can be used as an estima- 
tor of the total galactic mass, which in turns depends critically 
on the radial extent of the rotation curve. Therefore, after du- 
tifully noting that the regression analysis produces consistent 
fits whether the entire sample, or just the galaxies for which 
R(Vc)/R25 > 1 are used (Table 2), we will only retain the latter 
sample in the remainder of this paper. 

Because of its apparently negligible intrinsic scatter (a re- 
duced of 0.70, Table 2) and the fact that it does not follow di- 
rectly from simple dynamical arguments (see §3.1), the Vc-ctc 
relation is of obvious interest. However, by targeting only spi- 
ral galaxies it lacks generality; furthermore, it samples a narrow 
range in dc, ^ 90 to ^ 180 km s"'. It is interesting to compare 
the Vc-cTc relation observed for the spirals with the one de- 
rived recently by Kronawitter et al. (2000) and Gerhard et al. 
(2001) for a sample of 21 elliptical galaxies. For this sample, 
circular velocity curves are derived by applying non-parametric 
spherical models to the observed stellar absorption line profiles, 
velocity dispersion and surface brightness profiles. Kronawit- 

'NGC4486B is the only galaxy for which the models are deemed unreliable. 



ter et al. conclude that the rotation curves resulting from the 
models are very likely nearly self-similar: all rotation curves 
are consistent with being flat outside 0.3re, with the galaxy's 
effective radius^ As pointed out by Gerhard et al. (2001), 
is expected to be linearly related to any other velocity scale, for 
instance <Tc, as a consequence of the near dynamical homology. 
The constant of proportionality is likely to depend slightly on 
the details of the models; Kronawitter et al. (2000) and Ger- 
hard et al. (2001) do not discuss this point exphcitely, but by 
comparing their results to the those obtained by van der Marel 
(1991) and Lauer (1985) for the same galaxies, variations of the 
order of a few tens of a percent are possible if different dynam- 
ical models are used. 

Fig. 3 shows the Vc - ctc relation for the combined samples of 
16 spirals with R(vc)/R25 > 1 (soHd circles) plus the 20 ellipti- 
cals from Kronawitter et al. (2000, open circles). Values of a 
for the ellipticals are taken from Davies et al. (1987), and cor- 
rected to a common aperture size as detailed for the bulges. As 
expected, we recover the correlation between Vc and already 
reported by Gerhard et al. (2001). 

The fact that the spiral and elhptical samples follow the same 
relation is unexpected (see §3.1) and points to a common cause 
connecting galaxy dynamics at small and large scales across the 
entire Hubble sequence. A regression analysis to the ellipticals 
produces a fit consistent with the one obtained for the spirals 
once the two galaxies with ac < 70 km s~' are excluded (fits 3 
and 4 in Table 2). If a linear fit is forced, 0-^ = (0.65 ±0.02) for 
the ellipticals and (Tc = (0.61 ±0.01)Vc for the spirals. However, 
because the correlation between Vi,,> and CTc for the ellipticals is 
partly due to the scale-free nature of the models, we use only 
the sample of 13 spirals with /?(Ve)/^25 > 1 and CTc > 70 km 
s"^ in deriving the best fit to the Vc - <Tc relation: 

log = (0.84 ±0.09) log 0-^ + (0.55 ±0.19) (2) 

While this fit is strictly derived in the 70 < CTc < 180 km 
s"' range, it can also be considered valid in the 120 < CTc < 350 
km s~' range populated by the eUipticals (Fig. 3 and Table 2). 
As is the case for the M, - relation, the reduced Xr = 0-38 of 
the fit points to a relation with negUgible intrinsic scatter. 

3.1. Is the Vc — <^c Relation a Tautology? 

The tightness of the ~ o"e relation is suprising to the point 
of arising suspicions that we might be dealing with a tautology. 

The most trivial case to consider is the following: a strict pro- 
portionality between <Tc and Vc would result if both responded 
to the same mass distribution. For instance, if the region within 
which (Tc is measured extends to radii where the rotation curve 
is flat, ffc would naturally be linearly related to Vc simply by 
virtue of the virial theorem. This is not the case for any of the 
galaxies considered here: even for the bulges with the largest 
effective radius, ac only samples the innermost 0.5 kpc, while 
the rotation curve does not settle to a constant velocity until 
much farther out, typically 4 to 20 kpc. Fig. 4 plots the ra- 
tio of bulge central velocity dispersion to disk circular velocity 
against the ratio of the radii at which the two are measured (the 
radius at which Vc peaks is taken as a lower limit to the radius 
at which the circular velocity first settles onto a flat curve). As 
can be seen, R{Vmca)lR{(^c) ranges from values of a few to over 
100. 

Even so, if rotation curves of spiral galaxies formed a nearly 
homologous family, as is the case for the elliptical galaxies of 

and will not be considered further. 
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Fig. 4. — The ratio between bulge velocity dispersion (Jc and disk circular 
velocity Vc plotted against the ratio of the radii at which the two are measured. 
SoUd circles are galaxies for which i?(vc)/R25 > l O- 



rmendy & lUingworth 1983), a correlation between and CTc 
ensues. However, as in the case for the M, - (Tc relation, which 
is "expected" given the M, -Lb and the Faber- Jackson relation, 
the Ve-CTe relation is tighter, and therefore more fundamental, 
than any of the correlations mentioned above. 

One last possibility to explore is whether the - CTc relation 
could be nothing but the Tully-Fisher relation in disguise. Ver- 
heijen (2001) finds negligible intrinsic scatter for the /T-band 
Tully-Fisher relation when the circular velocity in the flat part 
of the rotation curve is substituted to the maximum rotational 
velocity. The tightness of the relation implies a fundamen- 
tal connection between DM halo mass and the total baryonic 
mass. The connection with bulge mass (and hence CTc) is how- 
ever not immediate (Norman, Sellwood & Hasan 1996). For in- 
stance, the bulge to disk fraction could conceivably depend on 
the detailed form of the halo angular momentum profile (Bul- 
lock et al. 2001) and the details of the angular momentum ex- 
change between gas clouds during dissipation (van den Bosch 
et al. 2002). It should be further noticed that Verheijen (2001) 
finds that the Tully-Fisher relation holds with negligible intrin- 
sic scatter down to Vc ~ 80 km s"\ while the Vc-CTc relation 
shown in Figure 3 seems to break down below Vc ^ 150 km 



(however, the constant of proportionality would depend on the 
details of the mass density profile, and there would be no rea- 
son for the elliptical and spiral samples to define a common 
relation). Unlike Gerhard's ellipticals, however, spiral galaxies 
show a strong dynamical non-homology: both the profile shape 
and the amplitude of the rotation curves depend on the galaxy 
luminosity (Casertano & van Gorkom 1991; Persic, Salucci & 
Stel 1996). Therefore, a hnear relation, or a relation at all, be- 
tween Vc and fTc is not expected a priori. 

The two arguments presented above also explain why the 
Vc-CTc relation is not a direct result of the so called "disk- 
halo conspiracy". Most rotation curves are characterized by 
a reasonably flat rotation curve in the outer parts. This impUes 
that the mass density profiles of the luminous and dark matter 
must be coupled, thereby "conspiring" to produce a flat rotation 
curve and the absence of conspicuous features marking the edge 
of the luminous component (Casertano & van Gorkom 1991). 
However the conspiracy breaks down in the inner, bulge dom- 
inated region, which is always characterized by a steeply in- 
creasing rotation curve whose detailed shape and overall slope 
depends strongly on the galaxy luminosity and morphological 
type (Persic, Salucci & Stel 1996; Sofue et al. 1999). For the 
galaxies in our sample, only for NGC 2841 the case could be 
made of a bulge-disk conspiracy: the rotation curve is approxi- 
matively flat all the way to the innermost point, ^ 2 kpc (Bege- 
man 1987). In all other cases, the bulge contribution is insignif- 
icant by the time the flat part of the rotation curve is reached 
(Broeils 1992; Kent 1987). In other words there is no direct ev- 
idence in these galaxies of a coupling between bulge and halo. 
At most, a connection between the Vc - (Tc relation and the disk- 
halo conspiracy can be thought of in the following terms. A 
rough correlation exists between maximum rotational velocity 
and Hubble type (e.g. Casertano & van Gorkom 1991); coupled 
with the correlation between Hubble type and bulge luminosity 
and between the latter and the bulge velocity dispersion (Ko- 

2nA=o.7, n„ = o.3,/j=o.7 



4. A RELATION BETWEEN SBHS AND DM HALOS 

Not being the result of a dynamical tautology, the near in- 
variance of CTc/ over two orders of magnitude in R( Vc) /R(cTc) 
seems to indicate a remarkable uniformity in the (luminous + 
dark) mass density profile along the Hubble sequence. In par- 
ticular, characterizing the slope and normalization of the rela- 
tion might help in constraining theoretical models and numer- 
ical simulation (e.g. Steinmetz & Muller 1995) following the 
formation and evolution of galaxies. 

Furthermore, as discussed in §2, the existence of the Vc - ctc 
relation is strongly suggestive of a tight correlation between 
M, and the total gravitational mass of the SBH host galaxy. 
If the rotation curves for all galaxies plotted in Fig. 3 are flat 
at a distance R from the center, i.e. v(R) = Vc, the virial theo- 
rem ensures that M(R) oc v/. Equations (1) and (2) then imply 
M, oc M(Kf-\ 

Several proposed scenarios trace the origin of SBHs to the 
very early stages of structure formation (Cattaneo, Haehnelt 
& Rees 1999; Silk & Rees 1998; Haehnelt & Rees 1993; 
Umemura, Loeb & Turner 1993). By defining the depth of 
the potential well, the mass and distribution of the DM halo 
control not only the formation of SBHs, but also their relation- 
ship to the luminous matter. Relating Vc to M^m is no trivial 
task, but at least a rough attempt can be made as follows. In a 
CDM-dominated universe. Mom is uniquely determined by the 
halo velocity v,,,> measured at the virial radius, R^r- The latter 
is defined as the radius at which the mean density exceeds the 
mean universal density by a constant factor, generally referred 
to as the "virial overdensity" Av,r. Based on this definition, 
and in virtue of the virial theorem, it immediately follows that 
Mdm must be proportional to the third power of Vi,,>. Being a 
function of Av,r, the constant of proportionaUty depends on the 
adopted cosmology and may vary with time (e.g. Bryan & Nor- 
man 1998; Navarro & Steinmetz 2000). For the purpose of this 
paper, we will adopt the same ACDM cosmological model used 
by Bullock et al. (2001)^. Under these assumptions: 
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Mdm = 2.7 X 10i2(Vv,>/200 km s-^f Mg 



(3) 




10» 10= 10' 10= 
SBH Mass M. (M^) 



Fig. 5 . — Same as Fig. 3, but with Cc transformed into SBH mass using equa- 
tion (1), and Vc transformed into DM halo mass following the prescriptions of 
Bullock et al. (2001). Solid circles are spirals, solid triangles ellipticals. Sym- 
bols accentuated by a larger open contour identify galaxies having a dynamical 
estimate of M,, which was used in the plot (Merritt & Ferrarese 2002). The 
upper limit on the SBH mass in NGC 598 is marked by the arrow. The dotted 
line represents a constant ratio Mom/M, = lO"". The solid line corresponds to 
the best fit using Bullock's et al. prescription to relate Vv,> to Vc (equation 5). 
The dashed line shows where the galaxies would lie if v^ir = Vc (equation 4), 
while if Vc/ Vvir = 1.8, as proposed by Seljak (2002), the points would move to 
the location shown by the dot-dashed line (equation 7). 



(Bullock et al. 2001). 

In a zeroth order approximation, we might assume that Vyj 
Vc. Combining equations (1), (2) and (3) gives: 



108 MfT, 



0.027 



DM 



12 



1.82 



(4) 



.0 vIQI^Mq 

However, for any reasonable shape of the DM density pro- 
file, the circular velocity will decrease towards 7?,,,>, in other 
words, Vc as listed in Table 1 is most certainly an upper limit to 
Vvir. The shape of the rotation curve, and therefore the Vc/v^r 
ratio is of course strongly dependent on the assumed DM den- 
sity profile. We will again follow Bullock et al. (2001) in 
adopting an NFW profile (Navarro, Frenk & White 1995) char- 
acterized by an "inner" radius r,. In this specific application, 
"Vvir I VciR) is a function of two variables, the concentration pa- 
rameter Cvir = Rvir/fs and the ratio R/rs (Bullock et al. 2001). 
N-body simulations (Bullock et al. 2001) show that for present 
day halos in the lO'^ - lO''' Mq range, Cv,> depends weakly 
on Mom- C^ir ~ 5 x \(f-M~f^^^, where Mv,r is in units of Mq. 
Therefore, for ~ a few x 10^ kpc (Klypin et al. 2001), 
rs = Rvir/Cvir ^ a few to a few x 10 kpc, i.e. is likely sam- 
pled by the flat HI rotation curves used in deriving the Vc - ctc 
relation. This is a reasonable conclusion: in the region of space 
dominated by DM, the turnover in the rotation curve can be 
roughly expected to coincide with the turnover in the DM den- 
sity profile. Therefore, by assuming R = r,, Vy,, / vdR) is a func- 
tion of Mum only. By using equation (3) to express Vv,,> as a 

'This is likely a lower limit to Vyir/ Vc(R), since the peak circular velocity is > 



function of Mom, Bullock's et al. results can be approximated 
as: 



Ml 



DM 



1012 



1.40 



200 km s"^ 



3.32 



(5) 



By using equation (5) to transform Vc to Mom, and equation 
(1) to transform tJc to M,, equation (2) can be written as (Fig. 
5): 



108 Mr 



0.10 



M, 



DM 



1012 Mr 



1.65 



(6) 



A third, completely independent way of converting the 
Vc — (Jc relation into a relation between SBH and DM halo 
masses is the following. Observational constraints on the mass 
density profiles of DM halos have been recently derived (Sel- 
jak 2002) by combining observations of galaxy-galaxy lensing 
with the Tully-Fisher and fundamental plane relations. For 
galaxies, the circular velocity at the virial radius decreases by a 
factor 1 .8 from its peak value at the optical radius. This decline 
is sleeper than predicted by Bullock's et al. (2001) numerical 
simulations, possibly due to the neglect, in the latter work, of 
the baryonic contribution to the rotation curves at optical radii 
(Seljak 2002). 

By assuming vdR) ~ 0.56^, the relation between Mom 
and M, becomes: 



M. 



108 Mr 



0.67 



© 



M, 



DM 



1012 M, 



© 



1.82 



(7) 



Comparing equations (4), (6) and (7) can give a rough idea of 
how the functional dependence of Mom on M, depends on the 
assumptions made for the DM density profile and the resulting 
rotation curve. 

4. 1 . The Significance and Scatter of the M, -Mom Relation 

The M,-MoM relation implies that the formation of SBHs 
is controlled, perhaps indirectly, by the properties of the DM 
halos in which they reside. Whether the connection between 
SBHs and DM halos is of a more fundamental nature than the 
one between SBHs and bulges, reflected in the M, - CTc relation, 
is unclear at this time. The M, - (7^ relation implies a scatter of 
the order of 20- 25% in M. for ac = 70-300 km s"i . The scatter 
in the M, -Mom relation is difficult to determine, since the un- 
certainty in the transformation between Vc and Mom is largely 
unknown. BuUock et al. (2001) estimate Alog(C,„>) = 0.18. In 
the range = 50-400 km s"', this translates into a 5% to 7% 
error in Mom in equation (5), and a 9% to 12% error on M, in 
equation (6). This indicates that, at least formally, the scatter 
in the M, -Mom relation could be a factor two smaller than the 
scatter in the M, - ac relation for any given choice of DM den- 
sity profile and cosmological parameters. Secure conclusions, 
however, will have to await an empirical characterization of the 
M, -Mom relation, with both M, and Mom determined directly 
from observations. 

Additional clues on how bulges and DM halos influence the 
formation of SBHs can be derived by studying bulgeless spi- 
rals. Such galaxies have been forcefully neglected in the present 
study, where we required the existence of a bulge velocity dis- 
persion in order to derive M,. Of the galaxies in Table 1, 
NGC598 is the one with the least conspicuous bulge"*. IC342 
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and NGC 6503 are late type spirals with minimal bulges; like 
NGC598, they deviate significantly from the relation defined 
by the more massive systems. All of the galaxies with no or 
minimal bulge for which a rotation curve exists, including both 
dwarf systems and low surface brightness (LSB) galaxies, have 
circular velocities below 150 km s~'(e.g. Blais-Oulette, Am- 
ram & Carignan 2001; de Blok & Bosma 2002). In particular, 
in the LSB compilation of van den Bosch et al. (2001) the only 
two galaxies with maximum circular velocity in excess of 150 
km s"' are also the only galaxies with a significant bulge com- 
ponent. Therefore, it does seem that bulges, as well as SBHs, 
are always formed in halos with circular velocity in excess of 
^ 150 km s"^ In less massive halos, the interplay between 
halo, bulge and the central SBH remains to be investigated. For 
this purpose attempts to directly measure the mass of a putative 
SBH in bulgeless galaxies is of primary importance. 

Whether equation (4), (6) or (7) is adopted, M, does not 
depend linearly on Mdm- According to equation (6), the ra- 
tio between SBH and halo mass decreases from ~ 2 x 10"^ at 
Mdm = 10"* Mq to - 10'^ at Mdm = 10'^ Mq, i.e. less mas- 
sive halos are less efficient in forming SBHs. This tendency 
becomes more pronounced for halos with Mdm ^5x10^' Mq : 
as mentioned above, NGC 6503, NGC 3198 and NGC 598 lie 
above the relation defined by the more massive galaxies. It must 
be pointed out, however, that the M, - ctc relation is not defined 
for M, < 10^ Mq, indeed we have no direct evidence of the 
existence of nuclear black holes in this mass range. If such 
SBHs exist, the recent report of an upper limit of only a few 
10^ Mq for the BH in NGC 598 (Merritt, Ferrarese & Joseph 
2001; Gebhardt et al. 2001), indicates that they are not more 
massive than would be inferred by extrapolating the M, - (Tc 
relation to the CTc ~ 20 km s~^ range. This excludes the possi- 
bihty that NGC 6503 and NGC 3198 could be reconciled with 
the M, -Mdm relation defined by the more massive galaxies. 

These findings are in qualitative agreement with the forma- 
tion scenario envisioned by Haehnelt, Natarajan & Rees (1998) 
and Silk & Rees (1998). In an attempt to match the DM halo 
mass function to the luminosity function of quasars, Haehnelt 
et al. noted that as the quasar lifetime is increased, a match can 
be obtained only at the expense of steepening the slope of the 
relation between M, and Mdm- In their model, a linear relation 
between M, and Mdm implies a quasar lifetime tgso ^ 3 x 10^ 
yr, about a factor 100 shorter than favored by current obser- 
vational constraints (e.g. Martini & Weinberg 2001; Ferrarese 
2002). If tQso ~ a few x 10^ yr, the slope of the M, -Mdm rela- 
tion must be increased to ^ 1 .7, in rather remarkable agreement 
with the results derived in this paper. 

Furthermore, Haehnelt et al. propose that SBH formation 
takes place during the very early stages of galaxy formation. 
Gas gravitating in a DM halo is deposited at the center of the 
halo at a rate which is proportional to Ihe Ihird power of the 
halo circular velocity Vy,>. Stars more massive than ~ 10* Mq 
undergo post-Newtonian instabiUties - and therefore collapse - 
before they can ignite hydrogen; for v„> larger than ~ 200 km 
s~\ such stars can be formed on short enough timescales that 
the collapse will take place before the star has had a chance to 
bum all of its fuel. Only halos with v^r ^ 200 km s"^ (or masses 
^ 10^2 Mq) would therefore host SBHs, and only SBHs with 



M, > 10*" Mq can be formed in such halos. It is therefore pos- 
sible that galaxies such as NGC 6503 and NGC 3198 do not 
contain SBHs. 

5. CONCLUSIONS 

The results in this paper can be summarized as follows: 

• We have considered a sample of 15 spiral galaxies with 
published measurements of both the bulge velocity dis- 
persion (Tc and the disk circular velocity Vc, measured 
beyond the optical radius. The sample displays a tight 
correlation between the two quantities. Furthermore, in 
a log((j)-log(Vc) plane, the narrow strip occupied by the 
spirals blends smoothly with that occupied by a sample 
of 20 elliptical galaxies (from Kronawitter et al. 2001 
and Gerhard et al. 2001). This is an interesting result, in 
view of the different methods used to derive for the two 
samples, and of the dissimilarities between the shapes of 
the rotation curves observed in spirals and modelled in 
elliptical galaxies. 

• The ellipticals and spirals define a common, linear (in a 

log-log plane) v,,- relation for 70 < 0-^ < 350 km s"' . 
Three galaxies have CTc below this range; all three deviate 
from the relation in the sense of having smaller than ex- 
pected values of CTc. The discrepancy is largest - a factor 
7 larger than the standard deviation on ac - for the galaxy 
with the smallest value of ac, NGC 598. 

• Rotation curves for all of the galaxies appear to be flat 
in the outer parts. In virtue of the M, - ac relation and 
the virial theorem, the Vc - ac relation must therefore im- 
ply that the masses of supermassive black holes correlate 
with the gravitational mass of their host galaxy, measured 
within a given radius R > R25. 

• If the further concession is made that, in CDM scenarios 
of structure formation, the virial mass of the DM halo is 
uniquely determined by Vc (Navarro & Steinmetz 2000, 
Bullock et al. 2001), the v^ — ctc relation is the observa- 
tional evidence that SBH masses correlate tightly with 
the mass of the DM halos in which, presumably, they 
are formed. The relation between M, and Mdm is non- 
linear: the ratio between M, and Mdm decreases with 
decreasing halo mass. This trend becomes more severe 
for Mdm ^ 5 x 10*' Mq. It is possible that less massive 
halos might indeed be unable to form SBHs, in qualita- 
tive agreement with the theoretical arguments proposed 
by Haehnelt, Natarajan & Rees (1998) and Silk & Rees 
(1998). 



I wish to thank Pat Cote, David Merritt, Milos Milosavljevic 
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'^The existence of a rarefied bulge in this galaxy is supported by both spectroscopic and photometric data (e.g. Minniti, Olszewki, & Rieke 1993; Regan & Vogel 
1994; Mighell & Rich 1995; Minniti 1996). 
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Table 1. Galaxy Sample^ 



Galaxy 


Hubble 


Distance 


Vc 


R{vc) 


R25 


(7 


0"c 


Ref2 


RefS 




Type 


(Mpc) 


(km s^^) 


(kpc) 


(kpc) 


(kill s^^) 


(km R^^) 


(Vr) 


{^) 








Galaxies with HI Rotation Gurves 








MW 


Sbl-II 


0.008 


180±20 


28.0 


11.4 


100±20 


100±20 


1 


1 


IC342 


SAB(rs)cd 


2.0 


185±10 


10.0 


14.0 


77±12 


77±12 


2 


2 


N224 


SA(s)b 


0.670 


240±20 


27.0 


21.4 


161±16 


146±15 


3 


3 


N598 


SA(s)cd 


0.85 


135±13 


15.5 


8.7 


27±7 


27±7 


4 


4 


N801 


Sc 


79.2 


216±9 


60.2 


37.3 


137±20 


144±27"' 


5 


5 


N2841 


SA(r)b 


18.0 


281±10 


81.2 


21.3 


192±13 


179±12 


6 


6 


N2903 


SB(s)d 


6.4 


180±4 


24.2 


11.7 


105±13 


106±13 


6 


7 


N2998 


SAB(rs)c 


67.4 


198±5 


46.7 


28.3 


107±24 


113±304 


5 


5 


N3198 


SB(rs)c 


9.4 


150±3 


46.1 


17.8 


65±10 


69±134 


6 


6 


N4258 


SAB(s)bc 


6.6 


210±20 


30.5 


19.5 


146±19 


138±18 


7 


1 


N4565 


SA(s)b 


17.1 


264±8 


51.2 


39.4 


151±13 


151±13 


8 


8 


N5033 


SA(s)c 


11.9 


195±5 


35.5 


18.5 


133±10 


122±9 


6 


7 


N5055 


SA(rs)bc 


8.0 


180±5 


44.0 


18.5 


115±7 


103±6 


9 


7 


N6503 


SA(s)cd 


5.9 


116±2 


11.4 


6.1 


45±7 


48±104 


6 


9 


N7331 


SA(s)b 


14.9 


239±5 


37.3 


46.0 


150±14 


139±14 


6 


7 






Galaxies with Optical Rotation Gi 


Lirves 








IC724 


Sa 


103.7 


302±3 


31 


35.3 


229±17 


243±18"' 


10 


3 


N753 


SAB(rs)bc 


90.0 


210±7 


21 


32.9 


129±18 


121±17 


11 


3 


N1353 


SA(rs)bc 


27.0 


191±6 


5 


13.3 


85±11 


92±12 


12 


3 


N1357 


SA(s)ab 


26.2 


259±5 


9 


10.7 


121±14 


121±14 


10 


3 


N1417 


SA(r)b 


54.2 


240±9 


14.5 


21.2 


148±18 


148±18 


13 


10 


N1620 


SAB(rs)bc 


60.7 


250±4 


20 


25.5 


117±19 


123±20"' 


12 


3 


N2639 


(R)SA(r)A: 


66.5 


318±4 


12 


17.6 


193±13 


195±13 


10 


3 


N2742 


SA(s)c 


24.2 


173±4 


9 


10.6 


68±30 


63±28 


11 


7 


N2775 


SA(r)ab 


24.3 


270±4 


10 


15.1 


170±14 


162±13 


10 


3 


N2815 


(R)SB(r)b 


40.0 


278±3 


18 


20.2 


199±20 


199±20 


12 


11 


N2844 


SA(r)a: 


26.4 


171±10 


6 


5.9 


114±12 


113±12 


10 


3 


N3067 


SAB(s)ab 


25.1 


148±1 


6 


8.9 


82±28 


84±29* 


12 


10 


N3145 


SB(rs)bc 


61.0 


261±3 


23 


27.4 


161±12 


166±12 


12 


3 


IN OZZO 


oJ\[T) DC 


46.5 


261±11 


21 


27.5 


157±16 


1 £J Q _1_ 1 


12 


3 


N3593 


SA(s)0/a 


12.4 


101±4 


2.1 


9.5 


56±15 


56±15 


10 


3 


N4062 


SA(s)c 


13.1 


154±13 


9 


7.7 


97±8 


90±7 


11 


12 


N4321 


SAB(s)bc 


15.1 


216±6 


10 


16.3 


83±12 


83±12 


11 


3 


N4378 


(R)SA(s)a 


43.1 


308±1 


9 


18.1 


206±19 


198±18 


10 


3 


N4448 


SB(r)ab 


16.8 


187±3 


6 


9.5 


179±21 


175±21 


12 


13 


N4698 


SA(s)ab 


20.0 


252±5 


8.2 


11.6 


134±19 


131±19 


10 


10 


N7217 


(R)SA(r)ab 


21.9 


241±4 


9 


12.4 


185±18 


171±17 


12 


13 


N7541 


SB(rs)bc 


51.0 


179±1 


12 


25.7 


104±14 


109±15'' 


11 


3 


N7606 


SA(s)b 


42.1 


240±4 


30 


32.9 


136±23 


124±2l4 


12 


10 



^Golumns are: 1) Galaxy name; 2) Revised Hubble type; 3) Galaxy distajice, from Ferrarese et al. 1999 or 
derived from the galaxy systemic velocity and Hq = 75 km s"'^ Mpc~^; 4) Gircular velocity measured at the 
outermost radius (col. 5) of the rotation curve, with references given in col. 9; 6) Galaxy diameter at the i? = 25 
mag arcsec"^ isophote (from the RG3, de Vaucouleurs et al. 1991); 7) bulge velocity dispersion, with references 
given in col. 10; 8) bulge velocity dispersion corrected to an aperture of size re/8. 

^References for Vc- 1. Oiling & Merrifield 1998 2. Newton 1980; 3. Kent 1989; 4. Gorbelli & Salucci 2000; 5. 
Broeils 1992; 6. Begeman 1987; 7. van Albada 1980; 8. Krumm & Salpeter 1979; 9. Bosma 1981; 10. Rubin et 
al. 1985; 11. Rubin, Ford & Thonnard 1980; 12. Rubin, et al. 1982; 13. Persic & Salucci 1995. 

^Reference for a: 1. Merritt, & Ferrarese 2000; 2. Terlevich et al. 1990; 3. Whitmore & Kirshner 1981; 
4. Merritt et al. 2001 5. Schechter 1983; 6. Bottema 1988; 7. Heraudean & Simien 1998; 8. Kormendy & 
lUingworth 1982; 9. Bottema 1989; 10. Heraudeau et al. 1999; 11. Bottema, R. 1992; 12. Whitmore et al. 1984; 
13. Whitmore et al. 1979. 

*No measurement of the bulge effective radius Re is available. The correction to the bulge velocity dispersion 
was calculated by assuming a ratio between re and R25 equal to the mean ratio calculated from the Baggett et 
al. (1998) for galaxies of the same Hubble type. 
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Table 2 

Results of the Linear Regression Fits: Ve = auc+P 



Sample^^^ 


N 


a 




X? 


1. All spirals 


37 


0.781±0.061 


0.70±0.13 


1.9 


2. Spirals, 7?(Vc)/7?25 > 1.0 only 


15 


0.683±0.060 


0.89±0.12 


0.70 


3. Spirals, R(Vc)/R25 > 1.0 and 


13 


0.839±0.089 


0.55±0.19 


0.38 


excluding N3 198 and N6503 










4. All Ellipticals 


20 


0.94±0.11 


0.31±0.26 


0.66 


5. Samples 3 and 4 combined 


33 


0.892±0.041 


0.44±0.09 


0.51 



Units for ac and Vc are km s 
^M33 is excluded from all fits. 



